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Abstract 

Wavepacket propagation calculations are reported for the interaction of a Rydberg hydrogen 
atom (n = 2 — 8) with Cu(lll) and Cu(lOO) surfaces (represented by a Chulkov potential), in 
comparison with a Jellium surface. Both copper surfaces have a projected band gap at the surface 
in the energy range degenerate with some or all of the Rydberg energies. The charge transfer 
of the Rydberg electron to the surface is found to be enhanced for n values at which there is a 
near-degeneracy between the Rydberg energy level and an image state or a surface state of the 
surface. The enhancement is facilitated by the strong overlap of the surface image-state orbital 
lying outside the surface and the orbital of the incoming Rydberg atom. These calculations point 
to the possibility of using Rydberg-surface collisions as a probe of surface electronic structure. 
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I. INTRODUCTION 


A Rydberg atom has one electron in a very weakly bound orbital, with a mean orbital 
radius that scales with where n is the principal quantum number. When such an atom 
approaches a solid surface, the state is strongly perturbed at long range by helds originating 
from the surface, due to the van der Waals interaction between the atom and its own 
image dipole in the surface, and also to localized charges at the surface Ultimately 

the perturbation may lead to charge transfer of the Rydberg electron to the surface ~ the 
surface ionization of the Rydberg atom - which typically occurs at a distance of around 4n^ 
in atomic units (e.g., for n = 20 the surface ionization distance is around 80 nm). In practice 
the ionization does not occur at a specihc distance but with an increasing probability as the 
atom approaches the surface (or equivalently as a varying ionization rate as a function of 
distance from the surface). 

An important question to ask is how much the ionization probability associated with 
the charge transfer process depends on the nature of the surface. The target surface can 
be physically structured (e.g. a rough surface, a stepped or vicinal surface, a surface with 
adsorbates) or electronically structured (e.g. band-gap structure, or thin dims or adlayers 
with embedded quantized states). It may also have localized charges as discussed above. For 
the charge transfer from a Rydberg atom to a conducting metal surface, the Rydberg energy 
is generally degenerate with the conduction band of the metal. Hence the charge transfer is 
resonant, and it would not be expected that the band structure of the metal would have a 
signihcant effect on this process. So for example gold and aluminum would be expected to 
show broadly similar behavior. 

However there are other situations where the surface does not have a true continuum of 
states at the appropriate energy for charge transfer. For example, a thin metal him behaves 
as a quantized ID well in the dimension perpendicular to the surface, and as a free-electron 
system (unquantized) in the dimensions parallel to the surface (thus forming a 2D band), 
as opposed to a Jellium surface which is effectively unquantized in all dimensions. Usman 
et al. studied such quantization effects on the charge transfer dynamics of ground state H" 
ions (from the negative ion to the metal surface) |1]. It was found that when the electronic 
energy of the H“ is close to that of a ID quantized well state, the rate of charge transfer is 
enhanced by the resonance, and the resultant ionization probability as a function of distance 
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was found to exhibit multiple peaks due to the many quantized thin film states that the 
atomic energy level crosses when the H~ approaches the thin him [1]. 

The observation of a resonance effect close to the ID quantized level can be explained by 
the preferential direction of charge transfer, which is along the z-axis on which the saddle 
point in the effective potential occurs, and so the charge transfer is most efficient when 
there is minimum transfer of momentum in the direction parallel to the surface, i.e., when 
the Rydberg energy is close to the ID quantized state. An alternative perspective is that 
states with momentum parallel to the surface have high angular momentum with regard to 
the position of the approaching atom, and therefore electron transfer into these states is 
hindered by conservation of angular momentum. A similar quantized-state resonance effect 
has been studied by Gauyacq and coworkers for the charge transfer of H~ ions at a Cu(lll) 
metal surface, but in the context of surface and image states embedded in the projected 
band-gap EHa. rather than quantized ID well states for the case of thin films. 

In this article the effects of electronically structured surfaces are investigated for the 
charge transfer of the Rydberg hydrogen atom, with principal quantum number n = 2 — 8, 
at the projected band gap surfaces Cu(lOO) and Cu(lll). The advantage of studying a 
Rydberg system compared to ground state H~ is that the resonance effects in charge transfer 
can be explored over a wide energy range by varying the principal quantum number of the 
Rydberg atom. Additionally high-energy image and surface states that lie in the energy 
range of experimentally accessible Rydberg hydrogen atoms can be probed. Such states are 
not easily probed by other methods. As is shown here, resonances occur at several values 
of n and at different values for the two surfaces. Wavepacket propagation calculations are 
presented that illustrate these resonant effects. 


II. THEORY 

A. Surface states and image states 

Surface states and image states are closely related, having energies within an energy 
gap of the band structure associated with the direction of the surface normal {z). Their 
electronic wavefunctions tend to zero towards the bulk metal as well as towards the vacuum 
(T —)■ 0, Ze —±C) 0 , where Ze is the electron position), and they are therefore mainly confined 
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to the metal-vacuum interface. Intrinsic surface states arise from the cleavage of the bulk 
metal (translational symmetry perpendicular to the surface is lost), and are localized mainly 
in the surface atomic layer with a wavefunction that decays exponentially towards the 

vacuum. Image states arise from the Coulomb-like attractive image potential for an electron 
outside the metal surface iV{ze) = — where is the image plane position) 

and from the surface barrier created by a gap of available bulk electronic states in the metal. 

Image states are localized mainly in the vacuum region of the interface (or more precisely, 
where the electron position is outside the image plane position rao The energies of the 
image states form a hydrogen-like Rydberg series (setting the usual nuclear charge term Z 
as 1/4): 

pis ^ ^ /I'l 

16 2(ni^g + a)2 ’ ^ ^ 

where riimg is the image-state index and a is the quantum defect parameter for a given 
surface (for the Cu(lll) and Cu(lOO) surfaces studied below, a is approximately 0.02 and 
0.24 respectively mm)- The corresponding ID wavefunctions have the form |T2j 


'^ninig(^e) ^e-Rnimg,^=o(^e/4), 


( 2 ) 


where Rm^ ^;=o(^e/4) is the normalized hydrogenic (s-wave) radial wavefunction. Thus, 


image states can extend far into the vacuum (< >n®„,g= 6(?7,img -|- a)^ao), and in the 

context of the Rydberg-surface interaction there can be signihcant overlap with the Rydberg 
electronic wavefunction at a long distance from the surface. Thus, the resonance effects in 
charge transfer via image states are expected to be particularly signihcant, especially when 
compared with the surface states which reside primarily inside the surface. For conducting 
metals such as gold, the image states are degenerate with the conduction band and the 
strong coupling of the image states to the degenerate continuum means that the energies 
are extremely broad and the resonant states short lived. Thus it is only for band-gap 
semiconductors (such as the Cu(lll) and Cu(lOO) surfaces considered here) where such 
states are sufficiently narrow in energy and long lived to observe resonant charge transfer. 

For an uncorrugated surface, the electron moves quasi-freely parallel to the surface (p), 
such that each surface state and image state forms a 2D continuum of states with energy 


pIS,SS/L 'I _ pIS,SS I 
(^ 11 ) - + 


2p ’ 


(3) 
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FIG. 1: (a) Energy of the electronic states in the model Cu(lll), Cu(lOO) and free electron Jellium 
surface, as a function of electron momentum parallel to the surface, A:||. The 3D bulk states are 
hatched and shaded in grey, the surface state (SS) and surface resonance (SR) are plotted as full 
black lines, image states (IS) and image resonance (IR) are plotted as full red lines, and the H 
atom Rydberg energies (l/2n^ a.u.) are shown as green horizontal lines, (b) is the same as (a) but 
expanded in the energy range of the Rydberg states of an external H atom. Rydberg levels close 
to a surface- or image-state or resonance are colored orange. The gamma point T at = 0 is also 
marked. The diagrams are based on those from Ref. m- 
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where the second term is the dispersion energy, with electron momentum parallel to the 
surface /cy, and an electron of effective reduced mass p ~ 1 (in atomic units) [H]. Figure 
shows the projected bulk band structure and the calculated energies of the surface and image 
states of the Cu(lll) and Cu(lOO) surfaces [I3]. Note that the full 3D bulk states are shaded 
in grey, and the projected band-gap is unshaded. Surface and image states that lie outside 
the energy gap, and are degenerate with the 3D bulk continuum, are commonly referred to 
as surface and image resonances (rather than ‘states’) and are broadened by the coupling to 
the continuum. 

B. Wavepacket propagation calculations 

The details of the implementation of the wavepacket propagation approach for the work 
presented here are given in Ref. [15]. In brief, the initial electronic wavefunction at atom- 
surface separations of Dq ~ 6n^ is found through matrix diagonalization using a Lagrange- 
Laguerre Discrete Variable Representation (DVR) for the radial coordinate [16] , and a Leg¬ 
endre DVR for the angular coordinate, 9 (the metal surface is assumed to be flat and isotropic 
in dimensions parallel to the surface plane, thus the system is cylindrically symmetric and 
rendered two-dimensional). The radial part of the initial electronic wavefunction is then 
projected onto a Coulomb Wave DVR na where the time propagation takes place. 

At each time step, the electronic wavefunction at a later time (or distance from the 
surface) is computed through the action of the time evolution operator, 

Tei(r;f + At) = f), (4) 

where iLei(i';t) is the electronic Hamiltonian which contains the kinetic energy operator, 
the centrifugal potential, and potential terms due the Coulombic interaction with the ion 
core, the applied electric held and the interaction of the electron with the surface. The time 
evolution propagator is approximated by a symmetric split operator [TH] and the details 
are given in Ref. [15]. In most of these calculations, the surface is moved towards the 
H atom at constant velocity and, unless otherwise stated, the velocity is taken to be v± = 
3 X 10“^ a.u. = 656m s“^ which is comparable to the mean velocity component perpendicular 
to the surface in experiments in our laboratory [191120]. 

In order to impose the correct outgoing boundary condition and to remove the spurious 
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reflections of the wavefunction from the edge of the numerical grid, a complex absorbing 
boundary is included close to the grid edge |2I] ■ To study the extent of the Rydberg-surface 
charge transfer, the total electron density on the numerical grid bound by the absorbing 
boundary is monitored as a function of atom-surface separation, and its derivative with 
respect to atom-surface separation corresponds to the ionization probability at the given 
distance. 

The advantage of using the CWDVR is its suitability for treating the Coulomb potential 
(in particular the singularity at the origin) and for describing the wavefunction in the metal, 
with a dense distribution of grid points near the origin and becoming regularly spaced at long 
range. This, together with favorable scaling of the wavepacket time-propagation (which only 
involves matrix multiplications, Eq. |^, allows calculations over a range of principal quantum 
numbers to become feasible, albeit still very computationally demanding. 


C. One-electron pseudo-potential 


Following the work of references [a [la El E2], the one-electron pseudo-potential of 
Chulkov et al. [m is used to model the Cu(lll) and Cu(lOO) surfaces in the direction 
normal to the surface {z), while the electron is allowed to move freely in the direction 
parallel to the surface (p). The potential has the analytical form: 


hee(^e) 

V,{z,) 

V2{z,) 

Vsiz,) 

V4{Ze) 


Vl{Ze) + V2{Ze) + V3(^e 
^10 + COS 

—A20 + A2 COs [/?(^ e )] 
R3exp[-a(2;e - 2 : 1 )] 

exp[-A(2:e-^im)]-l 

4(-2e 2ini) 


V4{Ze) 


Ze<0 
0 < Ze < Zi 

Zl ^ Z^ Z[jra 

Zira ^ Zg , 


(5) 


where is the bulk interlayer spacing, Aiq, Ai, A 2 and f3 are the independent model potential 
parameters, and the values for the Cu(lll) and Cu(lOO) surfaces are given in Table 
The remaining parameters in Eq. are determined by requiring the potential and its hrst 
derivative to be continuous: 


A20 — A2 — Aio — Al R3 — —A20 — ^ 

= g a = (6) 

A = 2a 2 ,,„ = 2i - pii (fg) . 
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Us (a.u.) 

Hio (eV) 

(eV) 

A 2 (eV) 

1 0 

Cu(lll) 

3.94 

-11.895 

5.14 

4.3279 

2.9416 

Cu(lOO) 

3.415 

-11.480 

6.10 

3.7820 

2.5390 


TABLE I: Values of the parameters used in the one-electron pseudo-potential given by Eq. for 
the Cu(lll) and Cu(lOO) surface [TT] . 


The bulk potential and the width and position of the energy gap is described by Vi in Eq. 
while the metal-vacuum interpolation potential and the energies of the surface and image 
states are described by V 2 and V 3 . The V 4 term describes the long range image-potential. 
The form of the Cu(lll) and Cu(lOO) pseudo-potential, and the energies and wavefunctions 
of the surface state and image states calculated from diagonalization of the ID potential are 
shown in Fig. For comparison the calculations for these surfaces are benchmarked against 
a Jellium surface [21] and the parameters used are those for aluminum. 

III. RESULTS 

A. Resonant charge transfer to surface states and image states 

The computation of the wavepacket propagation for the electronically structured Cu(lll) 
and Cu(lOO) surfaces is much more demanding than the Jellium surface. This is due to the 
periodic nature of the pseudo-potential inside the metal, which requires a large number 
of radial grid points, and the large components of the electron momentum parallel to the 
surface in the resonant charge transfer (see below), which requires a large angular momen¬ 
tum basis. Typically for the calculations carried out below, the numerical grid consists of 
~ 1000 X 200 (radial x angular) grid points. Even so, for the charge transfers that are 
restricted by the band-gap and unaided by the resonance effects of surface or image states 
(see below), it is found that high frequency oscillations are superimposed on the ionization 
probability profiles, and an increase of radial and angular points would alter the form of 
the prohles. However, it was found that for the results presented below, the calculations are 
sufficiently converged that the average ionization distance, and the range of atom-surface 
separations spanned by the ionization probability prohles, do not change signihcantly with 
a larger radial or angular basis, and so are acceptable for the comparison carried out below. 
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FIG. 2: Energies and wavefunctions of the bulk, surface and image states for the Cu(lll) and 
Cu(lOO) surface, calculated from the diagonalization of the ID pseudopotential given by Eq. 
with a Sine DVR basis [23]. The energies and wavefunctions are the eigenvalues and eigenvectors 
respectively obtained from the matrix diagonalization. SS, SR, IS, IR labels the surface state, 
surface resonance, image state and image resonance respectively. 

Furthermore, it is shown below that resonant charge transfer to band-gap embedded surface 
states and image states can be observed directly from visual snapshots of the wavepacket 
propagation, confirming that the different onset of ionization distances is not an artefact of 
the nnconverged calcnlations. 

The calculations presented here are for the H atom for which the levels of a given n have a 
degeneracy (arising from I and mi) of Under the surface potential (assnmed cylindrically 
symmetric) the /-degeneracy of the Rydberg n-manifold is broken, and the Rydberg manifold 
splits and forms states which are polarized with respect to the surface normal [25] akin to 
the Stark effect. If a held is present perpendicular to the surface (as typically occurs in 
experimental measurements [12] , then the energy levels are also split at long range from the 
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FIG. 3: Table of energies of the band gap, surface and image states or resonances for the Cu(lll) 
and Cu(lOO) surfaces at the gamma point F (A:|| = 0 ), compared to those of the Rydberg states of 
the hydrogen atom. SS, SR, IS, or IR labels the surface state, surface resonance, image state and 
image resonance respectively. 

surface by the normal linear Stark effect, and as the atom approaches the snrface from long 
range these states correlate adiabatically with the snrface-polarized states. To focns on the 
resonance effects of the band-gap embedded surface states and image states, the calculations 
are performed here on the most surface-oriented state (the most red-shifted state) of the 
n-manifold, which has the largest overlap with the surface. All the calculations reported 
here used an incoming velocity of 656 ms“^ for the Rydberg atom along the surface normal. 

The energies of the unperturbed H atom states are shown in Fig. [^as horizontal lines to 
compare with the energies of Jellium, Cu(lll) and Cu(lOO) surface electronic states. The 
values of the energies are given in Figure]^ to provide a more quantitative comparison and 
these suggest that there might be two or three resonances between the image or surface states 
and the H atom levels in the range n = 2 — 8. Figure shows a comparison of the mi = 0, 
n = 2 — 8 wavepacket propagation ionization probabilities calculated for a Jellium, Cu(lll) 
and Cu(lOO) model potential, as a function of scaled atom-surface separation (scaled by 
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Scaled atom-surface separation, D (n^ a„) 


FIG. 4: Ionization probabilities as a function of scaled atom-surface separation for the n = 2 — 8 
most surface-oriented Rydberg states calculated from the wavepacket propagation approach for a 
Jellium, Cu(lll) and Cu(lOO) model potential. Note that the n = 5 ionization probabilities are 
plotted twice to aid the comparison of the scaled ionization distances for each column. Rydberg 
states that are close to resonance with the surface or image states of Cu(lll) and Cu(lOO) surfaces 
are shown in bold. 


n“^). 

For n = 2, the charge transfer occurs at the greatest distance for Cu(lOO), then Jel¬ 
lium, and then for Cu(lll) at the shortest distance. Figure |^a) shows a snapshot of the 
wavepacket calculation for the three different surfaces; the atom is at a distance of 20.76ao 
from the surface. The amount of electron flux towards the surface clearly follows the order 
Cu(lOO) > Jellium > Cu (111). The n = 2 Rydberg energy level is nearly degenerate with 
the surface resonance of the Cu(lOO) surface at k\\ = 0 (F point, see Fig. 0 and Fig. |^, 
and so efficient charge transfer to the surface can occur via the surface resonance with little 
momentum transfer parallel to the surface - as illustrated by the directed beam of electron 
density moving towards the surface. The resonant effect is so strong that the charge transfer 
is even more efficient than for the free-electron Jellium surface. For the Cu(lll) surface, 
the n = 2 energy is in the energy gap at k\\ = 0, and so charge transfer to the surface state 
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FIG. 5: Snapshot of the electronic wavefunction, |'I'(r;t)p, during the wavepacket propagation 
calculation for the Jellium potential (left panels), Cu(lll) potential (middle panels) and Cu(lOO) 
potential (right panels) for the (a) n = 2, (b) n = 3 and (c) n = 4 Rydberg states. 

or the 3D bulk states can only occur with large k\\ component, which is inefficient since 
the dominant and preferential charge transfer axis is along z where the saddle point in the 
potential lies. Thus little charge transfer has occurred at this point of the trajectory - in 
accord with the ionization probability prohle shown in hgure 4. 

Figure shows that for n = 3, the charge transfer occurs much closer to the metal 
surface for both the copper surfaces compared with the Jellium. This is because the Rydberg 
energy level is in the energy gap of both Cu(lll) and Cu(lOO) surfaces at the T point (see 
Fig.|I] and Fig. |^, and so charge transfer occurs relatively inefficiently with non-negligible 
II components. This behavior is also apparent in hgure 5(b). For n = 4, the charge transfer 
distance is similar for the Cu(lll) and Jellium surfaces, while it is much smaller for Cu(lOO) 
surface. The snapshot of the wavepacket calculation in Figure |^c) for the three different 
surfaces illustrates that the more signihcant degree of ionization for the Jellium and Cu(lll) 
surfaces. The n = 4 Rydberg energy level is nearly degenerate with the image state of the 
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FIG. 6: Snapshot of the electronic wavefunction, during wavepacket propagation cal¬ 

culation for the Jellium potential (left panels), Cu(lll) potential (middle panels) and Cu(lOO) 
potential (right panels) for the (a) n = 5, (b) n = 7 and (c) n = 8 Rydberg states, n = 6 is not 
shown as the snapshot is qualitatively similar to n = 7. 

Cu(lll) surface, and so charge transfer can occur via the population of the image state, 
which can be clearly seen in the horizontal line of electron density, indicated by the arrow 
in the snapshot in Fig. [^c). For the Cu(lOO) surface, the Rydberg energy is within the 
band-gap and so charge transfer is inefficient. 

Above n = 4 the H atom energies are above the upper edge of the projected band-gap 
of the Cu(lll) surface, and are degenerate with the 3D bulk states (see Fig. 0 and Fig. 1^, 
and so the calculated ionization distance is similar to the Jellium model, as shown in Fig 
4 (although resonant transfer to image resonances can still be observed by studying the 
evolution of the electronic wavefunction, see below). For n = 5, the Rydberg energy is 
nearly degenerate with the hrst image state of the Cu(lOO) surface, and in the projected 
bandgap, (see Fig. and Fig. [^, and the population of this image state in the charge 
transfer process can be clearly seen in the snapshot of the wavepacket propagation shown 
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FIG. 7: The surface ionization probability of the n=4 H atom Rydberg state interacting with (a) 
Cu(lOO) and (b) Cu(lll) calcnlated using the constant velocity approximation. 

in Fig. g a). The resonance charge transfer to the image state results in a similar ionization 
distance to the Cu(lll) and Jellium model surface. This resonance effect to the hrst image 
state of Cu(lOO) is dost’ for n = 6 and 7, which involves a large k\\ transfer, and so the 
ionization distance becomes progressively smaller than the Cu(lll) and Jellium surfaces. 
For the Cu(lll) surface, ionization of the n = 7 Rydberg state via the close-in-energy image 
state resonance (see Fig. and Fig. g can be observed in Fig. gb). In fact, the effect of the 
image state resonance is so large that the charge transfer occurs at slightly larger distances 
than the Jellium surface (Fig. g. 

For the Cu(lOO) surface, the n = 8 Rydberg energy is close in energy with the second 
image state embedded in the projected band-gap (see Fig. and Fig. g. The population 
transferred to the second image state is illustrated in Fig. gc), and the resonance effect is 
reflected in the ionization probability shown in Fig. g (the ionization distance is similar to 
the Cu(lll) and Jellium model surface). 

B. Velocity effects 

The wavepacket calculations presented so far in this paper have used a constant velocity 
for the incident H atom. Experimentally we can alter the velocity of the incident H atoms 
by changing the delay between the pulsed excimer laser that produces H atoms from NH 3 
and the excitation laser pulses that produce the Rydberg atoms further down the supersonic 
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beam path [19]; the velocity can also be varied by altering the seeding gas. Using unseeded 
ammonia as the precursor gas we can access a range of perpendicular collisional velocities 
for the Rydberg H atoms from 500 m/s to 850 m/s while maintaining a usable molecular 
beam density. Figure shows repeats of the n = 4 WPP calculations for the Cu(lOO) and 
Cu(lll) surface at these typical experimental velocities. In general it would be expected 
that the higher the velocity, the less time the atom has to ionize in a given incremental range 
of positions D D — 6D. Therefore the probability distribution for the ionization distance 
should be shifted to shorter distances at higher velocities. Second if there are energy level 
crossings between the Rydberg atom and the surface image states that come into resonance 
over a specihc range of distances, then the crossing would be traversed more quickly at higher 
velocities and hence the process may be less sensitive to resonance at increasing velocity. 
The results shown in Fig. illustrate these expected trends. For the Cu(lll) surface in 
particular, there is a very clear shift of the probability distribution with velocity, with the 
peak probability shifting from d.ln^oo at 500 m/s to around S.Sn^ao at 850 m/s. The shift 
in the probability distribution is much less pronounced for the Cu(lOO) surface, as this is a 
non-resonant case for n = 4 H atoms, whereas it is resonant for Cu(lll). 

In the experiments described in [191 EO] there is an additional effect of the velocity im¬ 
pacting on the measurements associated with the probability that the ions produced are 
extracted and dectected after ionization. For greater collisional velocities the minimum dis¬ 
tance of ionization from which the ions can be extracted for a given held increases, such 
that the detection region shifts further from the surface. Coupled with the shift in ion¬ 
ization probability to smaller mean atom-surface separations, the integrated ion detection 
probability will decrease with increasing collisional velocity. 

In actuality, the velocity of the H atom towards the surface may not be constant, as 
there is an acceleration due to image-charge attraction. There are effectively three phases 
of acceleration: Initially on approaching the surface the Rydberg atom is accelerated by the 
van der Waals interaction of the neutral Rydberg atom to the surface, which scales as 
In the second phase of the wavepacket model, as the atom starts to ionize it develops an 
overall partial charge as the electron flux into the surface is generated, and hence there is 
a greater acceleration due to the image charge attraction to the surface. Finally the atom 
becomes fully ionized and there is a classical acceleration of the ion core due to its own image 
charge. These effects are neglected in the constant velocity calculations but can be included 
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FIG. 8: The surface ionization probability of the n=4 H atom Rydberg state interacting with (a) 
Cu(lOO) and (b) Cu(lll) calculated including the mean field term and in the presence of an applied 
field. 

by the addition of a mean field term [151 ESj. Here the forces acting upon the ion core are 
determined, and then the expectation value of the electronic Hamiltonian is calculated on the 
fly. The acceleration will tend to shift the probability of ionization towards lower distances. 
However, in order to model the experiments a further effect must be included; that of the 
applied held which is present for extracting the ions. This will oppose the acceleration of 
the partially charged atom in the second phase of acceleration. In these calculations we use 
very large (and experimentally impractical) helds, such as would be necessary to detect ions 
in experiments with n = A Rydberg atoms - in practice experiments are performed with 
much higher n values and hence much lower, more realistic helds are used. 

As shown in hgure inclusion of the mean held term, and the resulting acceleration, 
causes the expected reduction in the mean atom-surface separation for ionization. Addition¬ 
ally the diherence between ionization behavior for diherent incident velocities is reduced, 
because the acceleration tends to dominate the velocity ehects. For the on-resonance case, 
such as the image charge resonance at Cu(lll) for n = 4, the acceleration is less because the 
ionization occurs further from the surface (and hence image charge ehects are less) - thus, 
the ehects of the initial velocity are still more visible. 

The conclusion from the velocity-dependence calculations is that measurements of the 
velocity dependence of the surface ionization behavior may provide helpful additional indi¬ 
cation of resonant behavior. Recent experimental results from our group indicate that this 
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is indeed the case [20] , 


IV. CONCLUSION AND FUTURE PERSPECTIVE 

The surface ionization dynamics at an electronically structured surface are predicted by 
these calculations to be very different compared to a free-electron metal surface. Cu(lll) 
and Cu(lOO) have surface and image states embedded in a projected band-gap along the 
surface normal, forming a series of 2D bands. The electron motion inside the metal is 
restricted (based on the energetically available electronic states) in the direction normal to 
the surface, and is free in the direction parallel to the surface. Since the charge transfer 
of the Rydberg electron to the metal occurs preferentially along the surface normal (where 
the saddle point lies), with a minimum component of momentum parallel to the surface, 
the ionization distance is larger when the Rydberg electron is nearly degenerate in energy 
with one of the surface and image states. The resonant charge transfer to the surface and 
image states can be seen directly from the electronic wavefunctions during the wavepacket 
propagation. 

The charge transfer of Rydberg H atoms with a Cu(lOO) surface has been studied experi¬ 
mentally in our recent work [20], with principal quantum numbers in the range n = 25 — 34, 
but it has not been possible to study the charge transfer process at such low n values as 
reported in the calculations described here, owing to the short lifetimes of states and large 
helds required to extract and detect ions. On the other hand, extending the wavepacket cal¬ 
culations beyond n = 10 would be extremely challenging because of the very large basis sets 
required. Nevertheless it is possible to make some qualitative predictions for the higher-n 
experiments. 

For the Cu(lOO) surface, the projected band-gap extends to above the vacuum level 
(upper edge is at -1-2.98 eV), and the Rydberg series of image states, described by Eq. 
with a quantum defect of a = 0.24 (and a pre-factor of 1/15.93 instead of 1/16) [I2], will 
continue throughout this range. Thus, resonance effects are expected to play an important 
role in the charge transfer dynamics even for the experimental range of principal quantum 
numbers of n = 20 — 40. 

The image state energies of the Cu(lOO) surface in zero-held can be calculated from Eq. 
with the appropriate quantum defects However, it is well known from STM experiments 


17 



FIG. 9: The energies of the H atom Stark states and the energies of the Cu(lOO) image states as 
a function of electric field. The broad width of the Rydberg levels shows the approximate shift in 
energy from D = oo to D = 3n^ ao calculated using perturbation theory m- 

that the energies of image states can be strongly pertnrbed by external electric helds pTHSU] . 
Following Ref. |2B], the Stark shifted energies of the image states are calcnlated here from 
the diagonalization of the Chnlkov one-electron psendo-potential given by Eq. with the 
inclusion of a homogeneous electric held beyond the image plane. The calculated energies 
of the most red shifted n = 20 — 33 Stark states of the hydrogen atom, and the rii^g = 5 — 9 
image states of the Cu(lOO), as a function of the electric held are shown in Fig. |^b). The 
broad width of the Rydberg energies shown in Fig. j^b) represents the approximate shift 
in energy of the state due to the Rydberg-surface interaction from D = oo to D = 3n^ ao 
calculated using hrst order perturbation theory [31]. Figure [^b) shows that due the Stark 
shift of the Rydberg states, most Rydberg levels traverse at least one image-state resonance 
as the held increases, but that for n = 21, 25, 29,33 resonances are expected near zero held 
with image states IS5 to IS8. k = 0 Rydberg states, which do not shift in energy with 
external electric held, may provide the best experimental probe of these high lying image 
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states, as they provide the opportunity for the greatest overlap with the Copper localized 
image states. 

To our knowledge, to date, there has been no direct probe of the njm > 4 image states of 
Cu(lOO); in femtosecond two-photon photoemission experiments (where energies of < 3 
have been measured), the bandwidth of the laser excited a coherent superposition of high 
image states (riim > 4) to form wavepackets rather than exciting individual levels [32]. Thus 
the predictions shown above could lead to some exciting prospects for using Rydberg states 
as energy specihc probes of electronic states at surfaces. 
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